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ABSTRACT 

We use archival Hubble Space Telescope observations of resolved stellar populations 
to derive the star formation histories (SFHs) of eighteen nearby starburst dwarf galax- 
ies. In this first paper we present the observations, color-magnitude diagrams, and 
the SFHs of the eighteen starburst galaxies, based on a homogeneous approach to the 
data reduction, differential extinction, and treatment of photometric completeness. We 
adopt a star formation rate (SFR) threshold normalized to the average SFR of the 
individual system as a metric for classifying starbursts in SFHs derived from resolved 
stellar populations. This choice facilitates finding not only currently bursting galaxies 
but also "fossil" bursts increasing the sample size of starburst galaxies in the nearby 
(D< 8 Mpc) universe. Thirteen of the eighteen galaxies are experiencing ongoing bursts 
and five galaxies show fossil bursts. From our reconstructed SFHs, it is evident that 
the elevated SFRs of a burst are sustained for hundreds of Myr with variations on 
small timescales. A long > 100 Myr temporal baseline is thus fundamental to any 
starburst definition or identification method. The longer lived bursts rule out rapid 
"self-quenching" of starbursts on global scales. The bursting galaxies' gas consumption 
timescales are shorter than the Hubble time for all but one galaxy confirming the short- 
lived nature of starbursts based on fuel limitations. Additionally, we find the strength of 
the Ha emission usually correlates with the CMD based SFR during the last 4 — 10 Myr. 
However, in four cases, the Ha emission is significantly less than what is expected for 
models of starbursts; the discrepancy is due to the SFR changing on timescales of a few 
Myr. The inherently short timescale of the Ha emission limits identifying galaxies as 
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starbursts based on the current characteristics which may or may not be representative 
of the recent SFH of a galaxy. 

Subject headings: galaxies: starburst - galaxies: dwarf - galaxies: evolution - galax- 
ies: individual (Antlia, ESO154-023, SBS1415+437, UGC 4483, UGC 9128, NGC 625, 
NGC 784, NGC 1569, NGC 2366, NGC 4068, NGC 4163, NGC 4214, NGC 4449, 
NGC 5253, NGC 6456, NGC 6789, NGC 6822, IC 4662) 



1. A Perspective on the Starburst Phenomenon in Dwarf Galaxies 



A starburst is defined as a short-lived intense period of star formation (SF) that is unsustainable 
over the Hu bble time due to the lim ited s upply of gas within a galaxy. Starburst galaxies were first 
classified by lSearle &: Sargentl (119721 ) and lSearle et al.1 ()1973l ). based on the blue colors produced by 
the massive stars formed during the burst. In the local u niverse, starbursts create approximately 



1998 



10% of the radiant energy and 20% of the massive stars ([Heckman et al 
20041 ). At z = 1, starburst characteristics are found in 15% of galaxies (jO'Connell 



Brinchmann et al 



2005), presumably 



attributable to th e greater amounts of gas typically present in young galaxies and increased galactic 
interactions (e.g., 



Condon et al. 



1982 



Kennicutt et al.lll987l : iTelesco et alJll988h . 



The starburst 's impact on a galaxy and the surrounding intergalactic medium (IGM) is pri- 
marily due to the consumption of gas that fuels the burst and the feedback from massive stars 
formed in the burst. For example, one of the important mechanisms is the ionizing radiation pro- 
duced and the mass loss experienced during the lifetime of the massive stellar population and the 
resulting supernovae (SNe). The mass loss from massive stars alters the chemical composition of 
the interstellar medium (ISM) and future stellar populations. In addition, the flux of energetic 
photons escaping the galaxy changes the ionization state of the surrounding IGM. It is possible for 
starbursts to generate the energy and material to fuel galactic "super-winds" which can escap e from 
the low potential wells of dwarf galaxies enriching the metallicity of the IGM (|Skillmanlll997l ). This 
has been sug gested as an explanation for the low metallic ity measurements of many local dwarf 



galaxies (e.g.. ISpaans fc Norman 



1997 



Romano et al.ll2006l ). Yet, without understanding the SFH 



of a galaxy, it is difficult to constrain the impact of a starburst. 

Starbursts have been observed in gas rich galaxies such as disk galaxies and low mass dwarfs 



(e.g.. iHeckman et al.lll998l ). It is unclear whether most dwarf galaxies experience st arbursts during 



their evolution or whether only a subgroup of dwarf galaxies ex hibit bursts (e.g., iMarlowe et al 



19971 : iMaver et al. 



2001 



van Zee et al.ll2001bl : iPasetto et al.ll2003l ). It is possible that for the dwarf 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data Archive at 
the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. 
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galax ies with starbursts there may be a duty cycle to the bursts (jMarlowe et al,lll999l ; iLee et al 
20091 . and references therein) making this a repetitive mode of star formation in these galaxies. 



The ionizing radiation flux created during a starburst makes bursting galaxies easily studied at 
optical, ultraviolet (UV) , far infrared, and Ha rest wavelengths. Numerous studies over the past 
decade have used these different wavelen gth observations to study the properties of both nearby and 
distant actively bursting systems (e. 



19971 : iHeckman et al.ll2005l : lOev et al.ll2007l : ILee et al 



Colless et al 



1994 ; iBabul fc Ferguson! 11999 ; iMeurer et al 



20091 . among many others) . The goals of these 



studies have been many, including characterizing the different properties of a currently observable 
starburst and quantifying the feedback of starbursts to the different components of the ISM and 
IGM. 

Most studies have been linked to either instantaneous properties of starbursts (i.e., UV and 
Ha studies) or st atistical properties of galaxy samples. T he latter are typically b ased on broad- 



band colors fe.g.. lMadau et al.lll996T). spectral indices (e.g.. iKauffmann et al.ll2003l ) or full spectro 



photometric modeling (e.g., iBruzual &: Charlotl 120031 ) on integrated galaxy light. However, with 
resolved stellar populations, we can derive complete SFHs allowing us to study the evolution of 
starbursts over hundreds of Myr within individual galaxies. Among the m any relevant studies on 
resolved stellar populations in individual galaxies are those of NGC 1569 by Angeretti et al. (j2005l ) 
who estimate the SFH over the last 1 — 2 Gyr using synthetic color-magnitude diagrams (CMDs) 



and by iGrocholski et al.1 (|2008l ) who estimate the SF activity by overlaying metallici ty specific age 



isochrones on the obser ved CMP of NGC 15 69. A similar approach was taken by I Cannon et al 



(|2003h on NGC 625 and lAnnibali et~aTl (hoosh on N GC 4449 where SFRs were estimated by over- 
laying age isochrones. In the case of ICannon et al.l (|2003l ). the recent SFR estimates were refined 
in increments of 100 Myr using a model of blue helium burning (BHeB) stars. 

The high-resolution data used in such studies presents an opportunity to reconstruct the de- 
tailed SFHs over a longer temporal baseline for a sample of starburst dwarf galaxies in a homoge- 
neous and consistent way. To that end, we have gathered the observations of eighteen galaxies from 
the HST public archive and undertaken a detailed analysis of star formation on these systems. 
Using information from the CMDs and stellar evolutionary models, we derived SFRs for these 
galaxies as a function of time taking into account both systematic and statistical uncertainties. 
Our emphasis on uniformity allows us to consider the strength and importance of a starburst event 
in the context of the host galaxy's star formation in past epochs, and facilitates a direct comparison 
between galaxies while minimizing systematic uncertainties. 

In this first paper, Paper I, we present the observations and data reduction ($2]) and CMDs 
and SFHs ($3]) for eighteen dwarf starbursts. In §3] we discuss the classification of the galaxies as 
starbursts using a gas consumption timescale, compare our SFRs with existing measurements of 
Ha emission, and present a SFR normalized to the average SFR of each galaxy. The hypothesis 
that SF in a burst is "self-quenching" is discussed in £}5] and a summary of our conclusions in Sj6) In 
Paper II (jMcQuinn et alJIsubmittedf ) we will present additional analysis on these starbursts, and 
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derive the burst duration for each system. Our third paper in this series (McQuinn et al. in prep) 
will explore the spatial structure within the starbursts. 



2. Galaxy Observations and Photometric Reduction 
2.1. The Galaxy Sample 



Our sample of galaxies was selected from observations in the HST public archive and is listed 
in Table [TJ The majority of the galaxies were previously identified as starbursts in the literature 
using a number of methods. The different methods all rely on evidence of recent star formation such 
as the presence of high surface brightness regions, high Ha emission, or the blue color typical of 
massive, young stars (see Table [1] for citations and identification methods used). We explore a new 
method for identifying galaxies as starbursts based on the presence and distribution of young main 
sequence and helium burning stars in the optical CMDs. The higher concentration in dwarf galaxies 
of these young stars is a sign of recent SF activity that is likely more intense than the average SF 
experienced in low-mass systems. Using the characteristics of a CMD also has the potential to 
identify older or fossil bursts which may or may not be producing significant Ha emission at the 
current epoch. The identification of fossil starbursts potentially greatly increases the number of 
starbursts in the nearby universe that can be identified and studied. Using this method, we added 
four additional galaxies to our sample as possible starbursts (Antlia dwarf, ESO 154—023, NGC 784, 
and UGC 9128). The classification of these systems as starbursts is explored and discussed below 
in © 

The dwarf galaxies in the sample cover a range of brightness, morphologies (i.e., SBdm, 
IAB(s)m, pec, IBm), and spatial extent, facilitating an exploration of the starburst phenomenon in 
a range of dwarf galaxy properties. While the majority of the gala xies have high Galactic latitudes 
and low foreground extinction (i.e., < 0.2: ISchlegel et al.lll998l ). two of the systems, NGC 1569 
and NGC 6822, have high foreground extinction (Ar of 1.9 and 0.6 mag respectively). If we as- 
sume that the foreground extinction is uniform across an individual galaxy, the extinction can be 
easily taken into account. However, it is more l ikely that the foregrou nd extinction has variations 
and is not a uniform gray screen (see ^3.21 and iMcQuinn et alj 12009). We chose to include these 
well-studied systems in the sample with the caveat that the high extinction values increase the 
uncertainties on the SFRs. 

Our starburst sample was selected from observations in the HST archive based on three pho- 
tometric properties. First, we require both V and I band images of a galaxy in the archival 
observations. Second, we set a minimum photometric depth of the I band observations of ~2 mag 
below the ti p of the red g i ant branch; a requirement for a ccurately constraining the recent SFH 



of a galaxy ([Dolphin 



2002 



Dohm-Palmer &: Skillmanll2002l ). Third, the galaxies were required to 



lie close enough such that their stellar populations were resolved by HST imaging instruments. 
These properties were chosen to ensure robust reconstruction of the SFH. Seventeen of the eighteen 
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dwarf galaxies in our sample meet these criteria. We include an additional, galaxy to the sample, 
SBS 1415 + 437, even though the photometry stops somewhat short (~ 0.5 mag) of reaching our 
required photometric depth despite long integration times (P.I. Aloisi, Program I.D. 9361). While 
it does not strictly meet our photometric requirements, SBS 1415+437 is a well-studied starburst 
with an interesting recent SFH. We include it in our sample for comparison purposes and although 
the ancient SFH derived for this system is uncertain, the recent SFRs are still well constrained. 
Four other data sets were considered that were at least one magnitude above our photometric depth 
cut-off (UGCA 290 P.I. Schulte-Ladbeck, Program I.D. 8122; NGC 1705 P.I. Tosi, Program I.D. 
7506; NGC 3077 P.I. Seitzer, Program I.D. 8601; I Zw 18 P.I. Thuan, Program I.D. 9400). Testing 
showed that the recent SFHs on these galaxies were not properly constrained by the shallower 
photometry and are therefore not included in this work. 



The field of view of the HST observations are represented as footprints on DSS images (IBonnarel et al 



2000) in Figure [TJ Regions of active and quiescent star formation can be seen in each of the im- 
ages in these inhomogeneous systems. The galaxies span a large range both in angular size (from 
NGC 6822 at 15.5' to NGC 6789 at 1.3') and in distance (ranging from Antlia at 1.3 Mpc to 
NGC 784 at 5.9 Mpc with the outlier SBS 1415+437 at 13.6 Mpc). The HST imaging instruments' 
field of view covers a large fraction of most galaxies in the sample, with the exception of NGC 6822. 
The HST areal coverage thus includes a sufficient fraction of the stellar populations to profile the 
SFH although the SFRs derived for NGC 6822 will be lower limits, as noted below. The exposure 
times and basic data on each galaxy in the sample are listed in Table [TJ 



2.2. Photometric Processing 



Twelve of the eighteen galaxies were observed using the Advanced Camera for Surveys (ACS) 
Wide Field Channel (WFC) on HST. The observations for three of the galaxies were snapshots 
taken on a single HST orbit (CRSPLIT = 2) (i.e., NGC 4068, NGC 4163, and IC 4662), while 
the observations for the remaining objects were obtained across multiple HST orbits. Photometry 
was performed on the pipeline pro cessed, cosmic ray cleaned images using the ACS module of the 
DOLPHOT photometry package ([Dolphin! 12000) for the ACS images that were cosmic-ray split. 
The other six galaxies observed with the Wi de Field Planetary Ca mera 2 (WFPC2) instrument 
were reduced using the pipeline developed by Holtzman et al. (j20od ). Briefly, the WFPC2 images 
were pre-processed with the standard WFPC2 tasks mask, which flags bad pixels as saturated using 
data quality files and crmask which cleans the i mages. Photom etry was performed as part of the 
pipeline using the HSTphot photometry package ([Dolphin! 120001 ) in the PSF mode. Four galaxies, 
NGC 2366, NGC 4449, NGC 5253, and NGC 6822, were observed with multiple pointings. In these 
cases, the photometry measurements from each image were combined by geometrically modeling 
any overlapping regions in the images and selecting the photometry from only one of the images in 
the overlapping region. 



The definitions and derivations of the parameters from the two photometry methods, DOLPHOT 
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and HSTphot, used in the selecting high fidelity point sources are identical ensuring that our final 
photometry lists from the two HST instruments and photometric packages were homogeneously 
processed and filtered. The photometry output of DOLPHOT and HSTphot includes a charac- 
terization of the quality of each point source including the sharpness of the source, the degree of 
cr owding, and th e quality of the fit, etc; a detailed explanation of these parameters can be found 
in Dolphin (j200ol ) . We selected point sources from the photometry that were characterized as well- 
recovered stars with a minimum signal-to-noise ratio of 5. In addition to these criteria, the sum 
of the squared sharpness values for the two wavelength observations were required to be < 0.39 
for an individual point source. We rejected cuts based on crowding, because our tests showed this 
additional filtering preferentially eliminated bright point sources in areas of active star formation 
— one of the areas of interest for studying starbursts. Note that while crowding is typically not as 
significant in low surface brightness galaxies, in regions where crowding is high, the derived SFR 
may still be a lower limit. 

Artificial star tests were performed to quantify the completeness and systematic uncertainty 
of the photometry in each of the images. The artificial stars were evenly distributed across the 
field of view at the range of magnitudes and colors that bracketed the measured photometry. 
We present the percentage of stars recovered across each magnitude bin for one of the shallowest 
observations (NGC 4214) and for one of the deepest observations (Antlia dwarf galaxy) in our 
sample in Figure [2"P1 The artificial star lists were filtered according to the same parameters as the 
photometry lists. Point sources not meeting the criteria were retained as unrecovered stars. 



3. Star Formation Histories 

3.1. Connecting CMDs and Stellar Evolutionary Populations to SFRs 

CMDs created from our photometry are shown in Figure El with galaxies ordered from the 
faintest to the most luminous. Typical photometric errors per magnitude bin are shown. The 
different photometric depths across all the observations range from an absolute magnitude of Mj ~ 

2.5 for NGC 6822 to Mj 2 for UGC 6456 and NGC 4214. The observations stop short of the 

depth of the red clump, with the exception of Antlia, NGC 1569, NGC 2366, and NGC 6822. With 
the exception of SBS 1415 + 437, all photometry reaches ~2 mag below the tip of the red giant 
branch (TRGB ) in the I band, as needed to provide constraints on the ancient SFRs (t> 8 Gyr) 



(|Dolphinll2002l ). although our focus remains on the recent SFH. 

Multiple stellar evolutionary populations are identifiable by eye in the CMDs in Figure El 
including the main sequence (MS), the red giant branch (RGB), the red clump, the asymptotic giant 
branch (AGB), and the BHeB and red helium burning (RHeB) sequences. For didactic purposes, we 



2 Note that SBS1415 + 437 has a lower completeness limit than NGC 4214 despite long integration times, but is 
not representative of our galaxy sample. 
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present an enlargement of the CMD of NGC 4068 with the different stellar evolutionary populations 
overlaid in Figure [H While all of these evolutionary populations are used to derive the SFH, the 
helium burning (HeB) stars offer the clearest differentiation of starbursts, and are thus of the most 
interest in our study. 

S tars on the HeB sequences have masses greater than ^3 M and are burning helium in their 
cores (jBertelli et al.l [1994J; iDohm-Palmer fc Skillmanl |2002| ) . The ages of stars in this evolutionary 
sequence is between 5 — 1000 Myr during which time the stars will turn off from the MS and 
migrate across the CMD from the RHeB sequence to the BHeB sequence and back again. HeB 
stars older than ^600 Myr m erge into the red clump; the actual age is somewhat metallicity 
depe ndent (IBertelli et alJll994T). The l ocation of a HeB s tar on the CMD depe nds uniquely on age 



(i.e., IDohm-Palmer fc Skillmanl 120021 : IWeisz et al. 



2008 



McQuinn et al.1 l2009l ) in contrast to the 



MS and RGB where multiple aged stars with different metallicities can occupy the same space on 
a CMD. Because the HeB lifetimes are relatively short, the presence of HeB stars indicates that 
a galaxy has experienced recent star formation. While this is true for both the RHeB and BHeB 
sequences, the distribution o f the BHeB stars in the CMD is p articularly informative due the larger 
number of BHeB stars (e.g., IDohm-Palmer fc Skillmanl I2002T ) . 



For example, a galaxy with a constant, low SFR over the past 500 Myr will have a sparse but 
smooth BHeB star distribution with the number of BHeB increasing at fainter magnitudes (i.e., 
lower masses). Similarly, a galaxy with a constant, but significantly higher SFR over the same time 
period, will have a similar distribution in the BHeB stars although the overall number of BHeBs 
will be higher at all magnitudes (e.g., IC 4662 in Figure [3]). In contrast, for a galaxy with a variable 
SFR over the past 500 Myr, the BHeB distribution will be clumpy at different magnitudes such as in 
UGC 9128 (Figure |3|) where an over-density of BHeB stars is identified at intermediate magnitudes. 
This effect will sometimes be accompanied by a truncation of the MS at bright magnitudes if SF 
has recently stopped. The properties of the BHeB population in both the distribution as well as 
the total number of stars can differentiate a galaxy with constant, low-levels of recent SF, from a 
galaxy with bursting levels of recent SF. 



3.2. Methodology for Reconstructing SFHs 



SFH recovery program s are sophisticated tools that use stellar evolutionary theories to recon- 
struct past rates of SF (e.g.jTosi et aDll9jd; [T o lstoy & Sahalll996l : ballart et allll996l : iHoltzman et~al 
19991 ; iDolphinl |2002| ; lHarris &: Zaritskvl l200ll ) . Thes e techniques have been maturing over the past 
decade and are well-established in t he literature (see lSkillman Gallartll2002l ; ISkillman et al.ll2003l : 
Skillmanl 120051 ; iTolstov et al.l |2009|, and references therein) . The primary inputs in reconstructing 
the past SFR are the photometry and artificial star recovery fractions (i.e., quantitative observa- 
tional uncertainties and incompleteness) coupl ed with the ste llar evolutionary models. We employ 
the numerical CMD fitting program, MATCH (|Dolphinll2002l ) to construct a synthetic CMD based 
on the observed CMD varying the ages of the stellar populations to find the best fit to the observed 
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CMD, using the stellar evolutionary models of iMarigo &: Girardil ([20070 . We assume the Salpeter 
single sloped p ower law initial mass function with a spectral index of —1.35 from 0.1 — 100 M 



( Sarpeterlll955T) and a bin ary fraction of 35% with a flat secondary mass distribution. A recent paper 



by iGogarten et al.l (j2009i ) includes a detailed description of the parameters and method employed 



here on similar data. 

We constrained the metallicity, Z(t), to increase as a galaxy evolves in time. This physically 
motivated constraint guides the metallicity evolution in the absence of observational constraints that 
would have been available with deeper photometric depths. We did not constrain the metallicity 
function for three galaxies (Antlia, NGC 2366, NGC 6822 jf], as the deeper photometry reaches the 
bottom of the red clump (NGC 2366) or two magnitudes below the red clump (Antlia, NGC 6822) 
providing enough information at recent times (t<l Gyr) to constrain the metallicity. We present a 
comparison of the SFHs for both the unconstrained and the constrained metallicity solutions for the 
deepest photometric data, Antlia, in Figure [5j The two solutions show excellent agreement in both 
the ancient and the recent time bins. For shallower photometry, the metallicity constrained and 
unconstrained solutions also showed good agreement with small deviations in SFRs seen in only a 
few of the time bins. We chose to use the constrained metallicity solutions in these datasets as this 
produces a more realistic metallicity evolution in the galaxies avoiding large jumps in metallicity 



over short time periods ([Williams et al.ll2010l ). The galaxies that showed the largest deviations 
were the galaxies with higher extinction or significant stellar crowding (i.e., NGC 1569). As an 
example, we present two examples of the cumulative SFH with the metallicity constrained and 
unconstrained for two galaxy with typical photometric depth of our sample, ES0154— 023 and 
NGC 784, in Figure El As seen in the plots, the cumulative star formation is equivalent within the 
uncertainties. 

Photometric errors and differential extinction can broaden features in a CMD. These effects 
are explicitly accounted for in the CMD fitting programs employed in this study. The photometric 
uncertainties are quantified both with the uncertainties in the photometric measurements and with 
a measure of the photometric incompleteness provided by artificial star tests. Extinction is a 
free parameter fit by the SFH recovery program. While these factors must be considered in any 
SFH derivation, they are of smaller consequence in the deriving recent SFHs for dwarf starburst 
galaxies. Photometric uncertainties are smaller on the relatively bright stars found in sites of recent 
SF. Foreground extinction is lower as the galaxies in our study are predominately found at high 
Galactic latitudes and internal extinction is lower in the low metallicity environments of the sample. 
Nevertheless, to ensure a robust treatment of these parameters, we studied different regions of the 
galaxies with a range in photometric errors, and allowed for varying foreground and differential 
reddening when fitting the CMD. We found that the SFRs derived were not significantly affected 
by varying photometric errors and low extinction levels {Ay < 0.5). Higher levels of differential 



3 While the photometry of NGC 1569 is nearly as deep, the high levels of internal extinction dictated constraining 
the metallicity to increase with time. 
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extinction is more problematic affecting the SFRs derived for galaxies such as NGC 1569 and 
NGC 6822 with known high Galactic and/or internal extinction. For these systems, our reported 
u ncertainties in the deriv ed SFRs are higher. A complete discussion of these three factors is given 
in iMcQuinn et al.l ((2009). As an additional test of our results, we derived the SFH of one g alaxy, 
NGC 4163, using a different SFH recovery program, IAC— pop (|Aparicio fc Hidalgo! l2009l ). We 



found the SFHs to be equivalent within the uncertainties between the two methods. 

To illustrate the recovery of the SFH, in Figure [7] we present a synthetic CMD created by 
MATCH for ESO 154 - 023, alongside the observed CMD on the same scale. The synthetic CMD 
represents a typical model fit to the data in our sample. The different evolutionary populations 
are well-reproduced in the synthetic CMD, in both the broad and fine features. The quality of 
the fit between the observed and modeled CMDs can be quantified with an effective x 2 parameter 
( Dolphin! 120021 ) . The x 2 measures the likelihood of the SFH derived fro m the model C MD to be 
the true SFH of the observed galaxy given our inputs and models (i.e., iDolphinl |2002| ) . The x 2 s 
per degree of freedom for our analysis are presented in the last column of Table [2] and ranged from 
1.12 — 2.67. The three galaxies with a x 2 value above two suffer from high extinction (NGC 1569 
and NGC 4449) and high crowding (NGC 1569, NGC 5253, and NGC 4449). To match these 
features, MATCH fits a distance modulus and total extinction value (i.e., foreground plus internal 
extinction) to each galaxy. In Table [21 we compare these modeled values with distance modu li 
obtained from the literature and Galactic extinction values measured by ISchlegel et al.l (|1998l ). 
The distance moduli are well matched within a few t enths of a magnitude . The model extinction 
values show a wider spread when compared with the ISchlegel et al.l (119981 ) foreground extinction 
measurements. In many cases the modeled value is higher as is expected as the model estimates 
not only the foreground extinction but also the extinction internal to the galaxy. 

Crowding tends to be less in low surface brightness dwarf systems than in spiral and elliptical 
galaxies. However, regions of higher stellar density presented markedly lower completeness limits 
than the surrounding lower stellar density regions in four galaxies (IC 4662, NGC 1569, NGC 4449, 
NGC 5253). We accounted for this disparity by separating these galaxies into regions of higher and 
lower stellar density. In each region, we conducted the artificial star analysis thus more accurately 
describing the completeness across the changing condition of the galaxies. The SFHs were derived 
for each region separately and were summed creating a solution for the entire field of view. As a 
test, we performed the same analysis with two galaxies that did not exhibit the same extremes in 
stellar densities (NGC 4068 and NGC 4163) and found that the summed results were equivalent 
to the SFH derived from the entire field of view within the uncertainties. We therefore derived the 
SFHs from the global fields of view for the remainder of the sample which did not exhibit the same 
extremes in stellar density. 
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3.3. Star Formation Histories 

The lifetime SFHs (i.e., SFR(i)) are presented in Figure [8l No comparison with previously 
published SFHs are presented as none exist yet in the literature. The errors include both systematic 
and statistical uncertainties derived from Monte Carlo simulations. The average time resolution 
for deriving the SFH is <51og(t) ~ 0.3. However, the SFRs at large look-back times are not well 
constrained without photometry that reaches the oldest MS turn-off. Hence, we averaged the SFR 
results from the oldest time bins to create one coarse time bin of 6 — 14 Gyr; loosely constraining 
the ancient SFH. The uncertainties in this older time bin are relatively small due to averaging over 
a longer time period. By averaging over a long time interval, we traded temporal resolution for 
more secure SFRs thus providing a baseline for comparison with current SFRs and a context for 
the intensity of the recent SF. A finer time resolution is used in deriving the SFH in the most 
recent ~ 500 Myr leveraging the excellent temporal information from the HeB stars. The SFRs 
derived for NGC 6822 are lower limits due to the limited spatial coverage of the observations. As 
noted earlier in §2 . 1 1 the ancient SFH for SBS 1415+437 is not well constrained due to the shallow 
photometry. 

The uniform photometry and SFH recovery allow us to directly compare the results between 
galaxies. To first order, the SFHs share some common characteristics. Most show significant star 
formation in the oldest time bin indicative of galaxy assembly and the formation of stars in systems 
rich in gas. The galaxies also show elevated levels of recent star formation, as expected for systems 
previously identified as starbursts or fossil bursts. 

Deeper inspection reveals, however, that the details of the SFHs vary from galaxy to galaxy. 
While elevated levels of SF are sustained over a large interval of time (<5t>few 100 Myr), there are 
variations and inhomogeneities in the SFR profiles on small temporal scales (<5t~ 10 — 20 Myr). 
There does not appear to be a "typical" starburst SFR(t) profile. The most notable differences 
in the recent SFRs are more readily apparent in Figure [9] where we present an expansion of the 
last gigayear of the SFHs for each of the galaxies. A number of systems show a sustained SFR 
over a period of a few hundred Myr (e.g., ESO 154-023, NGC 4068, and UGC 4483) while others 
exhibit fluctuations in their recent SFR (e.g., UGC 9128, NGC 1569, SBS 1415 + 437, NGC 2366, 
NGC 625, and NGC 4214). One galaxy, NGC 6822 shows the beginning of a burst with elevated 
levels of SF in the most recent 50 Myr. 

One of the advantages of using resolved stellar populations is the ability to identify and study 
galaxies whose bursts have ended. Five SFHs show a decline in SFR in the most recent times (i.e., 
Antlia, UGC 9128, NGC 625, NGC 4163, and NGC 6789) after experiencing elevated levels of SF 
for a period of a few hundred Myr. These fossil burst galaxies are excellent systems in which to 
measure the complete duration of a burst event (see Paper II), rather than a lower limit. Fossil 
burst galaxies also increase the number of nearby systems in which one can study the starburst 
phenomenon. The resolved stellar populations required for any detailed SFH history limits such 
studies to galaxies within 5 — 8 Mpc. We encountered the limitations of studying farther galaxies 



even with observations of long integrations times (e.g., SBS 1415+437). Identifying fossil burst 
systems increases the sample of starburst galaxies in the nearby universe reducing our dependence 
on studying systems that push the limit of our observational capabilities. 



4. When Is Star Formation a Starburst? 

SF activity within a burst is complex and consequently, there are many definitions in the 
literature to identify and quantify bursts. Here, we consider three metrics to identify starbursting 
systems and discuss their applicability/limitations. 



4.1. Gas Consumption Timescales 



A fundamental characteristic of a starburst is that the high rate of SF during the burst is 
short-lived compared to the lifetime of the host galaxy. A limit for how long this high SFR may 
be sustained can be determined by calculating how long it would tak e to turn all th e raw material 
fueling the burst into stars. If this "gas consumption timescale" (T^c. lRobertslll963l ) is significantly 



less than the Hubble time, then the burst must be unsustainable and thus must be short lived on 
a cosmic time scale, increasing the chance that the elevated SFR is temporary. 

One can calculate r gas by dividing the mass of the gas (i.e., the raw material that fuels the 
starburst) by the peak SFR (M yr _1 ) during the starburst event. However, measuring an accurate 
and relevant gas mass is difficult for these low mass dwarf systems. The molecular phase of th e ISM 
is rarely detectable in dwarfs (e.g., Ilsrael et al.lll995l ; iTaylor et al.lll998l : iBarone et al.l 12000 ) . and 
thus is unlikely to contribute significa ntly to the gas supply, unless the uncertain co nversion factor 
from CO to H2 is extremely high (e.g.. lWilsonlll995l : llsraellll997l : lBolatto et al.ll2008l ). Atomic gas is 
thus likely to dominate the gaseous ISM. Indeed, the cu rrent best estimat es suggest that the ratio 
of molecular gas mass to atomic gas mass is 0.3 ± 0.05 (jLeroy et al.ll2005l . and references therein) 
indicating that the atomic gas mass dominates. While molecular gas is needed for star formation to 
proceed, the atomic gas is what ultimately fuels the formation of molecular gas, and is the relevant 
quantity to consider for the gas consumption times of low mass systems. Yet, measurements of the 
atomic gas mass can also be problematic. The atomic gas is frequently far more extended than the 
star forming disk reach ing up to several Hohnb erg radii at a density of 10 19 atoms cm -2 for dwarf 



irregular galaxies (e.g.. iHuchtmeier et al.lll98ll ). Thus, the gas mass available for SF is lower than 
the measured Hi mass in a galaxy making T gas an upper limit. 

Using the extended Hi mass as the best available proxy for the gas fueling SF in the galaxies, 
we present the resulting gas consumption timescales in the absence of stellar recycling in Table [3j 
We have calculated T gas from the peak SFR in the last Gyr for each galaxy (column 3) and atomic 
mass measurements taken from the literature (columns 4 and 5) after multiplyin g the Hi mass 



by a factor of 1.4 to account for the mass of helium and heavy elements (e.g., iRobertsI 11963 
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Kennicutt et al.l 11994 ). Of the sixteen galaxies for which we have found Hi measurements, fifteen 
show consumption timescales significantly less than a Hubble time (0.44 Gyr< T gas < 5.2 Gyr) 
indicating that the peak SFRs measured in these systems cannot be sustained for the entire history 
of the galaxy. Thus, while the gas consumption timescales are only upper limits, these fifteen 
galaxies still fit this simplistic starburst metric. One galaxy, NGC 6822, has a gas consumption 
timescale longer than a Hubble time {T gas = 26 Gyr); however, not only is the atomic gas know n 
to be unusually extended (> 600 arcmin 2 ) relative to the optical disk (jde Blok fc Walter! l2000l ) , 



but the optical observations cover only a fraction (~ 17 arcmin 2 ) of the star forming disk (see 
observational footprint in Figure [T]). 

Our simple calculation supports the gas consumption timescale as a useful first order metric 
for confirming that fifteen of the sixteen galaxies fit this starburst rubric. However, the r gas 
metric provides only a coarse evaluation of these systems as starbursts, due to the dwarfs' spatially 
extended Hi halo that may not be available for SF. An additional drawback of this simplistic gas 
consumption timescale metric is that is does not provide any additional insight into the nature of 
the burst. 



4.2. Starburst Thresholds 

Starbursts are defined by high levels of SF activity, yet the exact level of SF that constitutes 
"bursting" is subjective. Setting an absolute SFR as a threshold to identify bursts would bias any 
results toward larger galaxies with intrinsically higher SFRs. For an unbiased metric, one must 
consider the SFR relative to what an individual galaxy has experienced in the past. Previously, 
Scalo (jl986l ) and Kennicutt ( 1998 ) formulated the birthrate parameter (b) which compares the 



current SFR with the the average SFR over the lifetime of a galaxy (i.e., b= SFR / <SFR>). 
McQuinn et al. (2009) modified this parameter using the average SFR during the last 6 Gyr, instead 
of the average SFR over the lifetime of the galaxy (i.e., b recen t = SFR / <SFRo_6 Gyr >)• This 
calculation decouples any SF activity during the initial assembly period of galaxies from the SF 
during later epochs providing a more relevant baseline in the current epoch. Note that 6 Gyr is still 
a substantial fraction of the lifetime of a galaxy, and thus is a robust measure of the average SFR. 
In addition, the normalized SFRs are more securely anchored with SFRs from the more recent time 
bins (see §3.3p . 

Using a b recC nt > 2 to identify a burst and defining the beginning and end points of a burst at 



the time when b reC ent = 1 (jMcQuinn et al,ll2009l ) , the duration of the starburst events in our sample 
last over a few hundred Myr, with the shortest duration being 450 Myr and the longest being 1.3 
Gyr. The accuracy of these measurements depends in part on our inherent time resolution of the 
affected SFH. The temporal resolution over the last few 100 Myr is fine enough that we can discern 
and measure fluctuations in the SFR over short timescales (from <5t25 — 50 Myr in the last 100 Myr 
to <5t~ 125 Myr 500 Myr ago). In contrast, older bins have coarser inherent time resolution (from 
<5t~ 350 Myr in the 650 — 1000 Myr time bin), giving less accuracy on the interval burst duration. 
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Coarse time bins not only reduce the accuracy of the burst duration, but also reduce our sensitivity 
to bursts. With larger time bins, short duration increases are not detectable if their durations are 
much shorter than the width of the time bin in our derived SFH. We are therefore most sensitive 
to starbursts that have occurred in the last 500 Myr where are temporal resolution is the finest. 

The average SFR for recent times (b recen t = 1) are drawn in gray across the SFHs in Figure [9j 
the width of the gray line represents the u ncertainties in b re cent = 1. We identify starbursts if 



the birthrate parameter is greater than two (jKennicutt et al,ll2005l ). All galaxies show significant 



levels of recent star formation that are well above b recen t = 2 in multiple time bins. In the case of 
NGC 6822, the high levels of SF occur in the last 50 Myr, indicating that this burst is likely just 
beginning. The four additional galaxies we included in our sample based on the HeB distributions 
in their CMDs (Antlia, UGC 9128, ESO 154 - 023, and NGC 784) are confirmed as starbursts. All 
show brecent > 2 over multiple time bins. 

Because of the long temporal baseline available with CMD studies, our analysis also identifies 
fossil bursts whose current SFR is below our starburst threshold (i.e., Antlia, UGC 9128, NGC 625, 
NGC 4163, and NGC 6789). Four of the five fossil bursts (Antlia, UGC 9128, NGC 4163, and 
NGC 6789) do not show the typical burst characteristics of intense Ha, UV, or optical emission. 
For example, looking at a current optical image of the lowest mass galaxy in our sample, the Antlia 
dwarf galaxy, it is not at all obvious that it has undergone a period of significant star formation in 
its recent past. The absolute values of the SFR in Antlia are the lowest in our sample, but taken 
in the context of its own SFH, the recent burst of SF is significant and measurable. 

The fifth fossil burst, NGC 625, is unique among our sample. Although NGC 625 produces 



strong radio continuum and Ha emission often associated with starbursts ([ Cannon et all 120031 ; 



Cannon fc Skillmanl l2004h . the SFRs over the past 500 Myr fall below our starburst threshold 
of 

brecent — 2. Yet, NGC 625 shows strong evidence at other wavelengths of a starburst. This 
discrepancy can be understood if we look at the lifetime SFH in Figure El The SFH shows a low 
level of SF for the first 11 Gyr without the high levels of SF often associated with initial galaxy 
assembly, followed by multiple bursts in SF over the past 3 Gyr. Thi s "suppression" of SF at early 



times is similar to the SFH reported for the dwarf galaxy Leo A (jCole et al.1 120071 ). Our b reC ent 
metric assumes that the high levels of SF associated with early stages of SF activity in a galaxy 
occurs at t>6 Gyr; NGC 625 does not fit this assumption but we include it for comparison to the 
others in our sample. Note that suppressed initial SF over a significant period of time is also seen 
in two other systems in our sample, NGC 2366 and NGC 4214. 

FigurefTOlpresents the peak b reC ent values from each burst against the gas consumption timescale 
for fifteen of the systems. The observational field of view for one galaxy (NGC 6822) did not cover 
the majority of the optical disk so the gas consumption timescale for this galaxy yields an upper 
limit of ~ 26 Gyr; it not plotted as it is off the scale of this figure. There is little correlation between 
the gas consumption timescales and peak SFRs in the sample illustrating the limitations in using 
T gas for more than a simple diagnostic for classifying a galaxy as a starburst. The extended gas 
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disks of galaxies may not be available for SF within the optical disk rendering a gas consumption 
timescale a limited criterion with which to evaluate the properties of starbursts. 



4.3. Ha Emission from Starbursts 

The SF activity in a burst creates a short— lived population of high mass stars that emits 
copious amounts of radiation. This strong emission can be observed and used as a starburst 
indicator. Specifically, the high mass stars formed in a starburst create high levels of UV emission 
(M> 3 M Q ) and Ha emission (M> 20 M Q ) over the lifetime of these stars. Our time-resolved SFHs 
afford an opportunity to compare to the SFHs derived from their resolved stellar populations with 
the different timescales of these starburst indicators. We will explore the emission and timescales 
from the UV emission in a later paper; several of these galaxies are targets for a GAL EX legacy 



progr am (P.I. Skillman, Proposal number 60026). Here, we compare the Ha emission from lLee et al 



(|2009l ) with our SFHs. Since the Ha emission is expected to be emitted over a timescale < 5 Myr, 
we compare the Ha fluxes to the SFRs in our most recent time bin of 4 — 10 Myr. We do not 
include younger time bins because the evolution of the most massive stars with lifetimes < 4 Myr 
are not calculated in the stellar evolutionary models employed in the SFH program. 



Lee et al.l (|2009l ) identify a galaxy as a starburst if it has an Ha equivalent width (EW(Ha)) 



of gr eater than 100 A. They suggest that this corresponds to a birthrate value of 2 or 3 ( Kennicutt 



1998). This EW(Ha) starburst metric can be directly compared with our starburst indicator 
brecent , as shown in Figure [TTJ The horizontal dotted line represents the delineation of a starburst 
in EW(Ha) and the vertical dotted line represent our delineation of a starburst (b recC nt > 2). 
Galaxies presently bursting are plotted as blue points while galaxies with fossil bursts are plotted 
in thicker, red points. The measured EW(Ha) and b recen t values are listed in Table El 

The comparison of these two starburst indicators in Figure [11] highlights both the usefulness 
and the limitations of the techniques. Seven galaxies have EW(Ha) > 100 A classifying these 
systems as Ha starbursts. These systems also meet our starburst threshold of a b reC ent > 2 in our 
most closely matched time bin of 4 — 10 Myr. Eleven galaxies have a measured EW(Ha) below 
the 100 A threshold and were not classified as Ha starbursts. These eleven galaxies can be broken 
into two groups when considering our SFH results. In the first group, seven of these systems had 
a measured b re cent < 2 in the last 4 — 10 Myr suggesting that the starbursts in these galaxies 
are currently weak or have ended (red points). This is congruous with lower Ha emission. The 
second group, however, contains four galaxies (NGC 4068, ESO 154+023, NGC 784, NGC 4449) 
that show b r ecent > 2 in the last 4 — 10 Myr. In other words, these four galaxies show recent 
starburst characteristics from analysis of their stellar populations but do not meet the Ha emission 
starburst threshold. Note that our most recent time bin overlaps with, but is not coincident with, 
the timescale of Ha. The SFR in this final group of four galaxies must therefore be changing 
on timescales shorter than a few Myr to explain the discrepancy. Fluctuations of SFR on short 
timescales are also seen in the SFHs in Figure [9] at recent times where our temporal resolution is 
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the finest. 

A broader observation can be drawn from Figure [TT1 Analysis of stellar populations gives a 
longer timescale to consider SF activity within a galaxy. From this perspective, the starbursts 
are events lasting at least a few 100 Myr. If the starburst durations we measure are typical 
of dwarf galaxies, then short term fluctuations of SF activity are as fundamental as the longer 
term characteristics of the SFR. The inherently short timescale of the Ha emission is limited to 
identifying starbursts based only on the most current SF activity which, as we see in Figure [TT| 
may or may not be representative of the recent SFH of a galaxy. Ha emission will miss identifying 
an important number of starburst systems. 



5. Are Starbursts "Self-Quenching"? 



Feedback from massive stars in the form of UV ionizing radiation and SNe impact the sur- 
rounding gas cloud from which stars are being formed. The idea that this feedback can regulate 
or even quench in situ SF yia the heating and m echanical disruption of the gas has been put forth 



by previous authors (e.g., iNishi fc Tashird l2000i ) . Higher metallicity clouds are more effective at 



cooling so the quenchi ng effect would theoretically be more important in low-metallicity systems 
(bmukai fc Nishilll999h . 



The "sel f-quenching" proce ss has been extrapol a ted t o SF on galactic scales increasing with 
galaxy mass (|Kaviraj et al.ll2007l ). iFerguson fc Babull (|1998l ) posit that feedback from a burst limits 
the burst duration to short timescales of ~ 10 Myr. Numerous other studies h ave reported that star - 



bursts last < 10 Myr based on: (1) fitting starburst m odels to UV spectr a (jTremonti et al 



2001 



(2) studying stella r clusters within starburst galaxies ([Harris et al.l 120041) ; and (3) emission from 
Wolf-Rayet stars (jSchaerer. Contini. fc Kunthl Il999l ; iMas-Hesse &: Kunthl Il999i ). to name a few. 
These short durations would seem to corroborate the idea of "self-quenching" der ived by calcula- 



tions balancing the kinetic energy output of SNe with the binding energy of the gas (jThornlev et al 



2000j). If starburst durations are comparable to or shorter than the crossing time, they are self- 
extin guishing explosions destroyed by their energy output in a non-equilibrium fashion (IMeurer 

2oooh . 



In Paper II, we measure the durations of starbursts in the sample galaxies from our time- 
resolved SFHs presented in Figures [8]— |9l The durations range from 450 ± 50 to 640 ± 190 Myr in 
fourteen galaxies and up to 1300 ±290 Myr in four galaxies. These longer durations are comparable 
to or longer than the rotational period of the h ost galaxies and are i n agreement with other studie s 
measuring durations on the order of 100 Myr (jCalzetti et al.lll997l : iMeurerl l2000l : iLee et al.ll2009i ). 
The SFH show that starbursts events last a few hundred Myr, much longer than the lifetime of O 
stars. Rapid "self-quenching" does not apply to starbursts in our sample. 

The shorter timescales cited in the literature above likely result from measuring only part of 
a burst event, corresponding to the lifetime of individual star clusters formed. This "flickering" is 
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seen in the dips and variations in the SFHs on timescales of 10 Myr and suggests the small scale 
SF may be disrupted or quenched via stellar feedback. SF activity on these short timescales is 
associated with small spatial scales, and can be significantly affected by local conditions, making 
it a stochastic process. The non-equilibrium energy and mass transfer creating the small scale self- 
quenching may trigger SF in adjacent areas, creating self-propagating bursts of longer durations. 
We will investigate SF propagation within a burst in Paper III of this series using the spatial 
information of the stellar populations, regional SFHs, and stellar evolutionary time sequences. 



6. Conclusions 

We have used optical data obtained from the HST data archive to reconstruct the SFHs of 
eighteen nearby starburst dwarf galaxies. The SFHs show large variations from galaxy to galaxy, 
and, while the galaxies represent a class of objects undergoing significant recent star formation, 
there is no "idealized" SFR profile that characterizes the starburst phenomenon. Not only are the 
changes in the SFR during a burst varied but the peak SFRs during the bursts span nearly three 
orders of magnitude in value indicative of the range in physical conditions present in each of the 
galaxies studied. Thirteen of the galaxies show ongoing starbursts and five present fossil bursts. 
The SFRs fall on a continuum between bursting and "regular" with the distinctions between the 
two dependent on the individualized history of each host system and the complicated feedback 
mechanisms created during the bursts. With this in mind, we examined three metrics for identifying 
starbursts with the following results and conclusions: 

• The gas consumption timescale provides a coarse evaluation of whether a galaxy's SFR is 
sustainable over a cosmic timescale. For the sixteen galaxies with Hi mass measurements, fifteen 
systems can be classified as starbursts based on their gas consumption timescales being significantly 
shorter than the Hubble time. The gas consumption timescales are upper limits as the atomic gas 
mass used in the calculation measure the mass in the extended Hi halo; much of this gas may not 
be available for SF in the optical disk thus limiting the usefulness of such a metric. The remaining 
galaxy (NGC 6822) has an upper limit of the gas consumption timescales of 26 Gyr calculated from 
SFRs derived from a limited field of view of the optical disk. 

• In fourteen galaxies, the app lication of the EW(Ha) criterion for identifying bursting SF over 



the past ~5 Myr (jLee et al.ll2009l ) agrees with our measured b reC ent values in our closest matching 
time bin of 4 — 10 Myr. In four galaxies, the b recen t values from 4 — 10 Myr ago indicate bursting 
SF while the EW(Ha) measurements are below the starburst threshold indicating that while the 
bursts are long-lasting events, the SFR can change on timescales of only a few Myr. If the longer 
durations of a few hundred Myr seen in this galaxy sample are typical of starbursts, this longer 
lived phenomenon is best studied at wavelengths correlating with comparable emission timescales; 
the short timescale of Ha emission will miss identifying bursting galaxies experiencing short-lived 
fluctuations / dips in the SFR. 
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• We compared how the SFR changes over time within each galaxy normalized with an average 
SFR from the past 6 Gyr (b reC ent=SFR / <SFR> 0-6 Gyr); avoiding a comparison of absolute values 
of SFR biased towards larger systems. The recent SFRs are higher than twice the average SFR 
over the last 6 Gyr (b reC ent > 2) in seventeen galaxies confirming that these galaxies are starbursts 
by this stellar population metric. The exception is NGC 625 which presents a unique SFH. The 
high SFRs seen in most galaxies in our sample in the ancient time bins are not seen in NGC 625. 
The SFH shows quiescent levels of SF until ~ 3 Gyr ago, after which time SF proceeded at elevated 
levels until the current epoch. This galaxy shows evidence of a starburst at other wavelengths 
although it does fit our burst metric because of the high SFR baseline over the last 3 Gyr. 

• The starburst events seen in the SFHs last at least 450 Myr on timescales comparable to or 
longer than the rotational period of the host galaxies ruling out "self-quenching" of these bursts 
from stellar feedback mechanisms. 

7. Acknowledgments 

Support for this work was provided by NASA through grants AR-10945 and AR-11281 from 
the Space Telescope Science Institute, which is operated by Aura, Inc., under NASA contract 
NAS5-26555. E. D. S. is grateful for partial support from the University of Minnesota. J. J. D. 
was partially supported as a Wyckoff fellow. K. B. W. M. gratefully acknowledges Matthew, Cole, 
and Carling for their support. This research made use of NASA's Astrophysical Data System 
and the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics and 
Space Administration. Finally, we would like to acknowledge the usage of the HyperLeda database 
( http://leda.univ-lyonl.fr [ ). 

We would like to thank the anonymous referee for a prompt and very helpful report. 

Facilities: Hubble Space Telescope 

REFERENCES 

Aloisi, A., van der Marel, R. P., Mack, J., Leitherer, C, Sirianni, M., & Tosi, M. 2005, ApJ, 631, 
L45 

Angeretti, L., Tosi, M., Greggio, L., Sabbi, E., Aloisi, A., & Leitherer, C. 2005, AJ, 129, 2203 

Annibali, F., Aloisi, A., Mack, J., Tosi, M., van der Marel, R. P., Angeretti, L., Leitherer, C, & 
Sirianni, M. 2008, AJ, 135, 1900 

Aparicio, A., k Hidalgo, S. L. 2009, AJ, 138, 558 

Babul, A., & Ferguson, H. C. 1996, ApJ, 458, 100 



-18- 



Barnes, D. G., k de Blok, W. J. G. 2001, AJ, 122, 825 

Barone, L. T., Heithausen, A., Hiittemeister, S., Fritz, T., k Klein, U. 2000, MNRAS, 317, 649 

Begum, A., Chengalur, J. N., Karachentsev, I. D., Sharina, M. E., k Kaisin, S. S. 2008, MNRAS, 
386, 1667 

Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F., k Nasi, E. 1994, A&AS, 106, 275 
Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., k Blitz, L. 2008, ApJ, 686, 948 
Bonnarel, F., et al. 2000, A&AS, 143, 33 

Bottinelli, L., Gouguenheim, L., Fouque, P., k Paturel, G. 1990, A&AS, 82, 391 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, Kauffmann, G., Heckman, T., k 
Brinkmann, J. 2004, MNRAS, 351, 1151 

Bruzual, G., k Chariot, S. 2003, MNRAS, 344, 1000 

Calzetti, D., Meurer, G. R., Bohlin, R. C, Garnett, D. R., Kinney, A. L., Leitherer, C, k Storchi- 
Bergmann, T. 1997, AJ, 114, 1834 

Cannon, J. M., Dohm-Palmer, R. C, Skillman, E. D., Bomans, D. J., Cote, S., k Miller, B. W. 
2003, AJ, 126, 2806 

Cannon, J. M., k Skillman, E. D. 2004, ApJ, 610, 772 

Cole, A. A., et al. 2007, ApJ, 659, L17 

Colless, M., Schade, D., Broadhurst, T. J., k Ellis, R. S. 1994, MNRAS, 267, 1108 

Condon, J. J., Condon, M. A., Gisler, G., k Puschell, J. J. 1982, ApJ, 252, 102 al. 2009, ApJS, 
183, 67 

de Blok, W. J. G., k Walter, F. 2000, ApJ, 537, L95 

de Vaucouleurs, G., de Vaucouleurs, A., k Pence, W. 1974, ApJ, 194, L119 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, R. J., Paturel, G., k Fouque, P. 
1991, Third Reference Catalogue of Bright Galaxies, Springer, New York, NY (USA) 

Dohm-Palmer, R. C, k Skillman, E. D. 2002, A J, 123, 1433 

Dolphin, A. E. 2000, PASP, 112, 1383 

Dolphin, A. E., et al. 2001, MNRAS, 324, 249 

Dolphin, A. E., 2002, MNRAS, 332, 91 



-19- 



Drozdovsky, I. O., Schulte-Ladbeck, R. E., Hopp, U., Crone, M. M., & Greggio, L. 2001, ApJ, 551, 
L135 

Drozdovsky, I. O., Schulte-Ladbeck, R. E., Hopp, U., Greggio, L., & Crone, M. M. 2002, AJ, 124, 
811 

Ferguson, H. C, & Babul, A. 1998, MNRAS, 296, 585 
Gallagher, J. S., Ill, & Hunter, D. A. 1986, AJ, 92, 557 
Gallart, C, Aparicio, A., & Vilchez, J. M. 1996, AJ, 112, 1928 

Gieren, W., Pietrzyhski, G., Nalewajko, K., Soszyhski, I., Bresolin, F., Kudritzki, R.-P., Minniti, 
D., & Romanowsky, A. 2006, ApJ, 647, 1056 

Gil de Paz, A., Madore, B. F., & Pevunova, O. 2003, ApJS, 147, 29 

Gogarten, S. M., et al. 2009, ApJ, 691, 115 

Grocholski, A. J., et al. 2008, ApJ, 686, L79 

Harris, J., & Zaritsky, D. 2001, ApJS, 136, 25 

Harris, J., Calzetti, D., Gallagher, J. S., Smith, D. A. & Conselice, C. J., 2004, ApJ, 603, 503 

Heckman, T. M., Robert, C, Leitherer, C, Garnett, D. R., k van der Rydt, F. 1998, ApJ, 503, 
646 

Heckman, T. M., et al. 2005, ApJ, 619, L35 

Heydari-Malayeri, M., Melnick, J., k Martin, J.-M. 1990, A&A, 234, 99 
Hodge, P. W. 1980, ApJ, 241, 125 
Holtzman, J. A., et al. 1999, AJ, 118, 2262 

Holtzman, J. A., Afonso, C, k Dolphin, A. 2006, ApJS, 166, 534 

Huchra, J. P., Geller, M. J., Gallagher, J., Hunter, D., Hartmann, L., Fabbiano, G., k Aaronson, 
M. 1983, ApJ, 274, 125 

Huchtmeier, W. K., Seiradakis, J. H., k Materne, J. 1981, A&A, 102, 134 

Huchtmeier, W. K., Gopal, K., k Petrosian, A. 2005, VizieR Online Data Catalog, 343, 40887 

Hunter, D. A., Gallagher, J. S., k Rautenkranz, D. 1982, ApJS, 49, 53 

Israel, F. P., Tacconi, L. J., k Baas, F. 1995, A&A, 295, 599 

Israel, F. P. 1997, A&A, 328, 471 



-20- 



Karachentsev, I. D., et al. 2002, A&A, 383, 125 
Karachentsev, I. D., et al. 2003, A&A, 398, 467 
Karachentsev, I. D., et al. 2006, AJ, 131, 1361 
Kauffmann, G., et al. 2003, MNRAS, 341, 33 

Kaviraj, S., Kirkby, L. A., Silk, J., & Sarzi, M. 2007, MNRAS, 382, 960 

Kennicutt, R. C, Jr., Roettiger, K. A., Keel, W. C, van der Hulst, J. M., & Hummel, E. 1987, 
AJ, 93, 1011 

Kennicutt, R. C, Jr., Tamblyn, P., & Congdon, C. E. 1994, ApJ, 435, 22 

Kennicutt, R. C, Jr., 1998, ARA&A, 36, 189 

Kennicutt, R. C, Jr., & Skillman, E. D. 2001, AJ, 121, 1461 

Kennicutt, R. C, Jr., Lee, J. C, Funes, J. G., Sakai, S., &: Akiyama, S., 2005, ASSL Volume 
329, Starbursts from 30 Doradus to Lyman Break Galaxies, eds., De Gris, R. & Gonzalez 
Delgado, R. M., Springer (the Netherlands), 187 

Lee, J. C, et al. 2008, Astronomical Society of the Pacific Conference Series, 396, 151 

Lee, J. C, Kennicutt, R. C, Jose G. Funes, S. J., Sakai, S., & Akiyama, S. 2009, ApJ, 692, 1305 

Leroy, A., Bolatto, A. D., Simon, J. D., & Blitz, L. 2005, ApJ, 625, 763 

Loose, H.-H., & Thuan, T. X. 1986, ApJ, 309, 59 

Lopez-Sanchez, A. R., Koribalski, B., Esteban, C, Garcfa-Rojas, J. 2008, ASSP, eds. Koribalski, 
B.S.; Jerjen, H., Galaxies in the Local Volume, Springer (The Netherlands), p. 299 

Lopez-Sanchez, A. R., Koribalski, B., Esteban, C, Popping, A., van Eymeren J., & Hibbard, J. in 
press, Star Forming Galaxies Workshop Proceedings 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., Steidel, C. C, & Fruchter, A. 1996, 
MNRAS, 283, 1388 

Marigo, P., & Girardi, L. 2007, A&A, 469, 239 

Marlowe, A. T., Meurer, G. R., Heckman, T. M., & Schommer, R. 1997, ApJS, 112, 285 
Marlowe, A. T., Meurer, G. R., & Heckman, T. M. 1999, ApJ, 522, 183 
Mas-Hesse, J. M. & Kunth, D., 1999, A&A, 349, 765 

Mayer, L., Governato, F., Colpi, M., Moore, B., Quinn, T., Wadsley, J., Stadel, J., k, Lake, G. 
2001, ApJ, 559, 754 



- 21 - 



McQuinn, K. B. W., Skillman, E. D., Cannon, J. M., Dalcanton, J. J., Dolphin, A., Stark, D., & 
Weisz, D. 2009, ApJ, 695, 561 

McQuinn, K. B. W., et al., 2010, ApJ, submitted, Paper II 

Meurer, G. R., Heckman, T. M., Lehnert, M. D., Leitherer, C, & Lowenthal, J. 1997, AJ, 114, 54 
Meurer, G. R. 2000, Massive Stellar Clusters, 211, 81 

Mendez, B., Davis, M., Moustakas, J., Newman, J., Madore, B. F., & Freedman, W. L. 2002, AJ, 
124, 213 

Nishi, R., & Tashiro, M. 2000, ApJ, 537, 50 

O'Connell, R. W., 2005, ASSL Volume 329, Starbursts from 30 Doradus to Lyman Break Galaxies, 
eds. De Grijs, R. k Gonzalez Delgado, R. M., Springer (The Netherlands), 333 

Oey, M. S., et al. 2007, ApJ, 661, 801 

Omukai, K., k Nishi, R. 1999, ApJ, 518, 64 

Pasetto, S., Chiosi, C, k Carraro, G. 2003, A&A, 405, 931 

Paturel G., Petit C, Prugniel P., Theureau G., Rousseau J., Brouty M., Dubois P., Cambresy L., 
2003, A&A, 412, 45 

Pelupessy, F. I., van der Werf, P. P., k Icke, V. 2004, A&A, 422, 55 

Romano, D., Tosi, M. k Matteucci, F., 2006, MNRAS, 365, 759 

Reynolds, R. J. 1984, ApJ, 282, 191 

Roberts, M. S. 1963, ARA&A, 1, 149 

Sakai, S., Ferrarese, L., Kennicutt, R. C, Jr., k Saha, A. 2004, ApJ, 608, 42 
Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. M. 1986, Fundamentals of Cosmic Physics, Volume 11, 1 

Schaerer, D., Contini, T., k Kunth, D., 1999, A&A, 341, 399 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Searle, L., & Sargent, W. L. W. 1972, ApJ, 173, 25 

Searle, L., Sargent, W. L. W., & Bagnuolo, W. G. 1973, ApJ, 179, 427 

Skillman, E. D., Terlevich, R., Teuben, P. J., & van Woerden, H. 1988, A&A, 198, 33 



- 22 - 



Skillman, E. D. 1996, The Minnesota Lectures on Extragalactic Neutral Hydrogen, 106, 208 

Skillman, E. D. 1997, Revista Mexicana de Astronomia y Astrofisica Conference Series, 6, 36 

Skillman, E. D., & Gallart, C. 2002, Observed HR Diagrams and Stellar Evolution, 274, 535 

Skillman, E. D., Tolstoy, E., Cole, A. A., Dolphin, A. E., Saha, A., Gallagher, J. S., Dohm-Palmer, 
R. C, k Mateo, M. 2003, ApJ, 596, 253 

Skillman, E. D. 2005, New Astronomy Review, 49, 453 

Spaans, M. & Norman, C. A., 1997, ApJ, 483, 87 

Stinson, G. S., Dalcanton, J. J., Quinn, T., Kaufmann, T., & Wadsley, J. 2007, ApJ, 667, 170 

Swaters, R. A., & Balcells, M. 2002, A&A, 390, 863 

Taylor, C. L., Kobulnicky, H. A., & Skillman, E. D. 1998, AJ, 116, 2746 

Telesco, C. M., Wolstencroft, R. D., & Done, C. 1988, ApJ, 329, 174 

Thornley, M. D., Schreiber, N. M. F., Lutz, D., Genzel, R., Spoon, H. W. W. & Kunze, D., 2000, 
ApJ, 539, 641 

Thuan, T. X., & Seitzer, P. O. 1979, ApJ, 231, 327 

Thuan, T. X., & Martin, G. E. 1981, ApJ, 247, 823 

Tolstoy, E., & Saha, A. 1996, ApJ, 462, 672 

Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371 

Tosi, M., Greggio, L., & Focardi, P. 1989, Ap&SS, 156, 295 

Tremonti, C. A., Calzetti, D., Leitherer, C. & Heckman, T. M., 2001, ApJ, 555, 322 
Tully, R. B., et al. 2006, AJ, 132, 729 

van Eymeren, J., Marcelin, M., Koribalski, B., Dettmar, R.-J., Bomans, D. J., Gach, J.-L., & 
Balard, P. 2009, A&A, 493, 511 

van Zee, L. 2001, AJ, 121, 2003 

van Zee, L., Salzer, J. J., & Skillman, E. D. 2001, AJ, 122, 121 

Walter, F., Taylor, C. L., Hiittemeister, S., Scoville, N., k Mclntyre, V. 2001, AJ, 121, 727 

Walter, F., Brinks, E., de Blok, W. J. G., Bigiel, F., Kennicutt, R. C, Thornley, M. D., & Leroy, 
A. 2008, AJ, 136, 2563 



-23- 



Weisz, D. R., Skillman, E. D., Cannon, J. M., Dolphin, A. E., Kennicutt, R. C, Jr., Lee, J., Sz 
Walter, F. 2008, ApJ, 689, 160 

Whiting, A. B., Irwin, M. J., & Hau, G. K. T. 1997, AJ, 114, 996 

Williams, B. F., et al. 2009. larXiv:0911.4121l 

Williams, B. F., et al. 2010, ApJ, 709, 135 

Wilson, C. D. 1995, ApJ, 448, L97 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Note. — Cols. (3) Note that there arc multiple citations for many of these galaxies identifying them by different methods as starbursts. We show but one example citation per galaxy. 
(4) and (5): R.A. and Decl. in ,12000 coordinates. Col. (12) Diameter of the projected major axis of a galaxy at the isophotal level 25 mag arcsec - 2 in the B-band Paturcl ct al. (2003) 
HypcrLcda database (http://lcda.univ- lyonl.fr). Col. (13) Schlegcl ct U] i 19981 ). Col. (14) Absolute magnitudes uncorrected for Galactic extinction. 

Ide Vauc oulcurs et al] 119911 ) 

b Local Volume Legacy (LVL) project, Lee ct al. (2008) 

References. — (1) | Aloisi et al.N2005l) ; (2) ICannon et al.N2003D ; (3) iDolphin et all i2001|); (4) iDrozdovsky et al l )2001|); (5) iDrozdovskv et al.l )2002i); (6) IGieren et al.1 12006]); (7) 
IGrocholski et al.l i200a ); < 8) iKarachentsev et al.H2002f) : (9) IKarachentsev et abU2003i) : flO) lKarachentsev et alj 120061 ); (lp lMendez et al.H2002l ); (12) ISakai et al.U2004l ); (13) lTullv et al.l 
120061 ); 
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Table 2. Comparison of Distance, Extinction, and Fit Values 
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Note. — Cols. (2) Distance Modulus best fit by the CMD fitting program. The 
uncertainties are lower bounds as they include only statistical uncertainties and (3) 
Distance Modulus reported by various authors in literature. See Table[l]for references. 
(4) Foreground and internal extinction best fit by the CMD fitting program. The un- 
certainties are lower bounds as they include only st atistical uncertainties . (5) Galactic 
foreground extinction at \r = 650 nm reported bv lSchlegel et ah! il998T) . (6) Best \ 2 
value of CMD fitting program. 

a 'HSB' refers to the high surface brightness region. 

b 'LSB' refers to the low surface brightness region. 
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Table 3. Classification of Starbursts 
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References. — (1) Begum et al. 2008 (FIGGS); (2) Cannon et al. 2004; (3) de Blok 2000; (4) Walter et al. 2008 
(THINGS SURVEY); (6) van Eymeren et al. 2009; (7) van Zee 2001; (8) Whiting et al. 1997.; (9) Thuan & Seitzcr 1979; 
(10) Lopez-Sanchez et al. 2008; (11) Lopez-Sanchez et al. 2009; (12) Huchtmeier et al. 2005; (13) Bottinelli et al. 1990; 
(14) Swaters & Balcells (2002) 

Note. — The gas consumption timescales (r gas ) were calculated using peak SFR and the Hi gas mass adjusted by a 
factor of 1.4 to account for helium and heavy elements; they are upper l imits as the atomic gas extends up to several 
Holmberg radii farther than the optical disk for dwarf irregular galaxies (Huchtme ier et al]|l98lh . The uncertainties in 
Tgas were calculated by adding in quadrature the SFR uncertainties with an assumed uncertainty of 10% in the atomic 
gas mass measurements. The EW(H«) measurements were taken from lLee et al.l (120091 ). 

a The gas consumption timescale for NGC 6822 is an extreme upper limit. Not only is the atomic gas unusually extended 
relative to the optical disk, but the field of view of the optical observations used to derive the SFRs does not cover the 
main optical disk. 
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Fig. 1. — The footprints of the HST archival observation fields of view u sed to derive the SF Hs 
shown on DSS V band images: Antlia, UGC 9128, UGC 4483, NGC 4163 faonnarel et alJhoonh . 
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Fig. m— HST footprints continued: UGC 6456, NGC 6789, NGC 1569, NGC 4068 
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Fig. m— HST footprints continued: NGC 625, NGC 784, NGC 5253, NGC 6822 




Fig. [Q— HST footprints continued: NGC 4214, NGC 4449 
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V and I Absolute Magnitudes 

Fig. 2. — The completeness functions for the V band (solid lines) and I band (dotted lines) for the 
shallowest photometry in our study (top panel: NGC 4214) and the deepest photometry in our 
study (bottom panel: Antlia). These completeness limits bracket the range in photometry for all 
eighteen galaxies with the exception of SBS 1415 + 437 whose completeness is lower than NGC 4214 
as noted in the text ( §2.11) . 
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Anilia UGC9128 




V-l V-l 

Fig. 3.— Color magnitude diagrams of Antlia, UGC 9128, UGC 4483, NGC 4163 with average 
photometric uncertainties per magnitude bin. The contours show the density of points where they 
would otherwise saturate the plot. Contour levels are 7.5, 15, 30, 60, 120, and 240 thousand points 
mag -2 . The distances used to calculate the absolute magnitudes are noted in Table [TJ The TRGB 
break occurs at Mj = —4.0. Note the differences in photometric depth and the presence and 
distribution of BHeB stars. 
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Fig. El— CMDs continued: UGC 6456, NGC 6789, NGC 1569, NGC 4068 
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Fig. El— CMDs continued: SBS1415+437, IC 4662, ESO154- 023, NGC 2366 
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Fig. 4. — CMD of NGC 4068 with the evolutionary stages of the stellar populations labeled. The 
MS, BHeB, RHeB, RGB, AGB, and red clump evolutionary stages are all easily identified in the 
stellar populations. 
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Fig. 5. — The best SFH fit with metallicity constrained to increase in time for the Antlia dwarf 
galaxy is plotted as a solid line with the best SFH fit with unconstrained metallicity evolution 
overplotted in shaded gray. The solutions are in excellent agreement in both the ancient and recent 
time bins. 
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Fig. 6.— The cumulative percentage of SF in ESO 154-023 and NGC 784 found both with the 
metallicity constrained to increase in time (black lines) and the metallicity unconstrained in the 
solutions (dashed gray lines). In both galaxies, the cumulative SF are in very good agreement. 
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Fig. 7. — The left panel shows the Hess diagrams for the ESO 154-023 observations. The right 
panel is the best-fit synthetic Hess diagram to the observations representative of a typical fit for 
the sample galaxies. While there are slight discrepancies, each of the evolutionary populations are 
well-described in the synthetic CMD. 
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Fig. 8.— The SFHs of Antlia, UGC 9128, UGC 4483, NGC 4163 derived from the HST optical 
observations of the resolved stellar populations for the eighteen galaxy sample. The oldest time 
time (6 — 14 Gyr) is an average SFR that constrains the ancient SF giving context for the overall 
SFH yet these numbers are not used in our calculations. The uncertainties in this age bin are 
underestimated due to averaging over a longer time. The time resolution at more recent times is 
much finer and makes use of the unambiguous age dating of the HeB stars. Note the elevated levels 
of star formation in these most recent time bins. Expansions of the last 1 Gyr SFHs are presented 
in Figure [H The galaxies are ordered by Mb luminosity from faintest to brightest. 
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Fig. El— Lifetime SFHs continued: UGC 6456, NGC 6789, NGC 1569, NGC 4068 







-44- 



0.2 
0.15 

0.1 
0.05 

0.08 
0.06 
^ 0.04 

L 0.02 

© 

w 0.1 
0.05 

0.15 
0.1 
0.05 



SBS 1415+437 



H r- 



IC 4662 



H h 



ESO 154-023 



H h 



NGC 2366 



j i i_ 



J_ 



-\ 1 — ^ 



1 - i T l — (- 



-\ 1 h 



H h 



H 1 h 



I n l T | | 



_l I I L 



10 5 
Age (Gyr) 



J I I L 







Fig. El— Lifetime SFHs continued: SBS1415+437, IC 4662, ESO154- 023, NGC 2366. Note the 
ancient SFH for SBS 1415+437 is not well constrained due to the shallow photometry. 
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Fig. H— Lifetime SFHs continued: NGC 625, NGC 784, NGC 5253, NGC 6822. Note the SFRs 
derived for NGC 6822 are lower limits due to the limited spatial coverage of the observations. 
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Fig. El— Lifetime SFHs continued: NGC 4214, NGC 4449 
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Fig. 9.— The last 1 Gyr of the SFHs of Antlia, UGC 9128, UGC 4483, and NGC 4163 presented 
in Figure [8] are shown here in greater detail. The gray line represents the average SFR over the 
past 6 Gyr (b recen t = !)• The thickness of the gray line represents the uncertainties in the average 
which are affected by the photometric depth, extinction, and photometric crowding. Each profile 
shows elevated SFRs at recent times yet the SFRs over the past few hundred Myr are diverse. All 
galaxies exhibit the common characteristic of having SFRs well above the last 6 Gyr average yet 
there is no specific pattern of SFR that characterizes a burst nor is there an absolute level of SF 
that sets a threshold for bursting SF. The galaxies are ordered by Mb luminosity from faintest to 
brightest. 
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Fig. El— Recent SFHs continued: UGC 6456, NGC 6789, NGC 1569, NGC 4068 
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Fig. [91— Recent SFHs continued: SBS 1415+437, IC 4662, ESO154-023, NGC 2366 
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Fig. El— Recent SFHs continued: NGC 625, NGC 784, NGC 5253, NGC 6822 
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Fig. El— Recent SFHs continued: NGC 4214, NGC 4449 
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Fig. 10. — The gas consumption timescale is plotted as a function of the maximum value of b re cent 
during the most recent burst. Gas consumption timescales are based on Hi mass measurements, 
which are available for all but two of the galaxies in our sample. We consider the gas consump- 
tion times to be upper limits, as it is unlikely that the significant fraction of gas at large radii 



(Huchtmeier et al. 1981 



The uncertainties in r gas were cal- 



will be available for star formation, 
culated assuming a 10% uncertainty in the atomic gas masses. The plot indicates that the gas 
consumption timescale is not a reliable indicator of the strength of a starburst, given the lack of 
correlation with b reC ent- On the other hand, all of the consumption times are less than a Hubble 
time, indicating that the recent SF cannot be sustained indefinitely. 
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Fig. 11.— The EW(Ha) (A) (ILee et alJkoog ) is plotted against the 

brecent values from our time 

bin of 4 — 10 Myr; the most closely matche d with the Ha em ission timescale. The horizontal dotted 
line represents the 100 A threshold set bv lLee et al.l ( 20091 ) delineating a starburst galaxy and the 
vertical dotted line represents our limit of b re cent = 2 for the same purpose based on resolved stellar 
population analysis. Seven galaxies classified as Ha starbursts also show bursting levels of SF in 
our most recent time bin. Seven galaxies that do not meet the Ha starburst threshold, also show 
lower levels of SF in the past 4—10 Myr. Four galaxies show discrepancies between the two metrics 
indicating that the SF is changing on short timescales of a few Myr. These galaxies clearly show 
bursting SF from a stellar populations perspective but do not show starburst levels of Ha emission. 



